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Abstract. The formation of image-potential states at the interface between a 
graphene layer and a metal surface is studied by means of model calculations. An 
analytical one-dimensional model-potential for the combined system is constructed and 
used to calculate energies and wave functions of the image-potential states at the f- 
point as a function of the graphene-metal distance. It is demonstrated how the double 
series of image-potential states of free-standing graphene evolves into interfacial states 
that interact with both surfaces at intermediate distances, and finally into a single 
series of states resembling those of a clean metal surface covered by a monoatomic 
spacer layer. The model quantitatively reproduces experimental data available for 
graphene/Ir(lll) and graphene/Ru(0001), systems which strongly differ in interaction 
strength and therefore adsorption distance. Moreover, it provides a clear physical 
explanation for the different binding energies and lifetimes of the first (n = 1) image- 
potential state in the valley and hill areas of the strongly corrugated moire superlattice 
of graphene/Ru(0001). 


PACS numbers: 71.15.-m, 73.20.-r, 73.22.Pr, 78.47.J- 

Keywords'. image-potential states, graphene, metal surfaces, electronic structure 



Formation of Image-Potential States at the Graphene/Metal Interface 

1 . Introduction 


2 


The ability to fabricate freestanding graphene, a single atomic layer of graphite, has 
raised interesting qnestions with regard to its snrface electronic structnre. Dne to the 
fact that this system exhibits a mirror plane with two snrfaces, the eigenstates of this 
system mnst be either a symmetric or antisymmetric snperposition of the electronic 
states of each snrface. This applies also to image-potential states, a class of intrinsic 
snrface states that exist at all solid snrfaces due to the interaction of an electron in the 
vacuum with the polarizable surface. This interaction can be described by the classical 
image-potential which gives rise to a hydrogen-like Rydberg series of electronic surface 
states and resonances ES- Consequently, it has been predicted that freestanding 
graphene should possess a double Rydberg-like series of image-potential states of even 
and odd symmetry [^. However, there are no experiments so far on freestanding 
graphene that provide clear evidence for these two series. 

Image-potential states have been observed when graphene forms an interface on 
metal or semiconductor surfaces where it can be grown with remarkable high quality and 
long-range order [^[^. Experiments on graphene/SiC(0001) using scanning-tunneling 
spectroscopy (STS) and two-photon photoemission spectroscopy (2PPE) have 
revealed a splitting at least of the hrst (n = 1) image-potential state. It has been 
argued that these states are remnants of the first even and odd state of freestanding 
graphene even if its mirror symmetry is in fact broken due to the presence of the 
substrate. In contrast to these findings, only a single series of image-potential states 
has been observed on weakly interacting graphene/metal systems 11-14 . For other 


metals, like Ni, Pd, Rh, and Ru, the interaction between graphene and the metal 
can be much stronger. Together with the lattice mismatch, this typically leads to a 
pronounced buckling of the graphene layer forming a moire superstructure with large 
periodicity [^[^. Graphene/Ni(lll) represents an exception with an almost vanishing 
lattice mismatch. Despite the strong interaction, this results in the growth of a flat 
graphene layer at a distance of dg = 2.1 A 16 . For g/Ru(0001), on the other hand, the 


corrugation is particularly pronounced and the hexagonal moire superlattice is divided 
into strongly interacting valley {L) areas with a distance of dg = 2.2 A and weakly 
interacting hill (H) areas with dg = 3.7 A 


17 


18 


The observation of held-emission 
resonances with different energies in the L and H areas by STS raised a controversial 
debate about the assignment of the lowest members of the series of Stark-shifted image- 
potential states 19 -22 . Time and angle-resolved 2PPE experiments 23 support the 


general assignment of Borca et al. 19 . Based on the 2PPE results, it has been suggested 


that the series of image-potential states is slightly decoupled from the Ru substrate in 
the L areas, whereas the hrst (n = 1') image-potential state in the H areas has a 
substantial amplitude below the graphene hills. This explains the larger binding energy, 
shorter lifetime, and higher effective mass of this interfacial state which is closely related 
to the interlayer state of graphite and of other layered materials [5 |^,25 . 


Different model potentials have been developed for the description of image- 
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potential states of freestanding graphene and graphene/metal systems. Full ab initio 
approaches are not suited for this purpose because the accurate description of the 
long-range image force is too costly. Silkin et al. constructed a hybrid potential for 
freestanding graphene, which combines a potential, derived from a self-consistent local 
density (LDA) calculation for the description of the short-range properties and an image- 
potential tail for the proper description of the asymptotic long-range properties [^. 
De Andres et al. used an alternative approach that is based on the dielectric response of 
the graphene and requires only a minimal set of free parameters . For the description 
of graphene/metal systems, however, either the metal 26 , the graphene layer 27 


or both has been represented only by simple effective barriers characterized by 
empirical parameters that are fitted to the experimental results on the binding energies 
of the specific system. 

Here, we present an analytical one-dimensional model-potential that makes it 
possible to calculate wave functions and binding energies of image-potential states at the 
F-point for different graphene/metal systems that are characterized by specific graphene- 
metal distances. For this purpose we combine potentials for a realistic description of 
the projected metal band gap, the image-potential of the metal, and the potential of the 
graphene for an arbitrary graphene-metal distance. Furthermore, we include the doping 
of the graphene layer by considering the work function difference between graphene and 
the metal as well as corrections due to higher-order image-charges. All parameters of the 
model potential are determined from the properties of the bare metal surfaces and the 
freestanding graphene and no free parameters are used for the calculations at different 
graphene-metal distances. This distance has a strong impact on the properties of the 
image-potential states. We show how the double Rydberg-like series of freestanding 
graphene evolves to a single series when a flat graphene layer approaches the metal 
surface as realized for the weakly interacting g/Ir(lll) and the strongly interacting 
g/Ni(lll), which exhibit distinct different graphene-metal separations. The observation 
of two series can be explained for corrugated graphene layers with different local 
graphene-metal separations as demonstrated for g/Ru(0001). The wave functions of the 
two series of this system are found to show significant differences in their symmetry, 
binding energy and probability density close to the Ru(OOOl) surface. The latter 
determines the electronic coupling to the metal bulk and can be connected to the 
experimentally observed difference of the inelastic lifetimes of the corresponding states 
with the same quantum number in the L and H areas. 


2. Analytical Model-Potential 

For the description of the graphene/metal interface we have constructed a one¬ 
dimensional model-potential 

V(z) = V^(z) + Vg{z - dg) + V^{z) + 6V{z), (1) 

which is composed of a metal potential and a graphene potential Vg{z), where 

dg is the spatial separation of the graphene layer in the ^-direction with respect to the 
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position of the last metal surface atom a.t z = 0. V$(z) and 6V{z) are corrections 
which consider the work function difference between graphene and the metal and the 
influence of higher-order image-charges, respectively. We note that our potential does 
not contain any free parameters for the fitting to experimental binding energies of the 
combined system. 

For the modeling of the metal, we use the well-established one-dimensional 
analytical potential introduced by Chulkov et al. 28 . This potential describes the 


bulk by the two-band model of nearly free electrons and matches the asymptotic image- 
potential in order to achieve a total potential \fn(z) that is continuously differentiable 
for all z. Its free parameters [see Appendix A| are fitted to experimental data on the 
size and position of the projected band gap as well as on the binding energies of the 
image-potential states of the clean metal. The limitations of the two-band model for 
the description of the d-bands in metals like Ru, Ir, and Ni are not important in the 
framework of the present work because the properties of the image-potential states are 
most sensitively dependent on the surface-projected bulk band structure in the vicinity 
of the vacuum level i^vac- 

The graphene layer is modeled by a parameterized analytical potential Vg{z) 
which has been fitted to the numerically calculated ”LDA-|-image tail” hybrid potential 
developed for the description of freestanding graphene by Silkin et al. [^. The binding 
energies of the image-potential states for freestanding graphene, which we obtain by 
using our analytical potential, agree well with those obtained with the use of the 
”LDA-|-image tail” hybrid potential [see Appendix B| . Although it is known that the a 
and TT bands of graphene are considerably shifted for strongly interacting graphene/metal 
systems 29,30 , we neglect this change of the electronic structure by using the same 
potential for all graphene-metal distances since this change has only a minor influence 
on the image-potential states at the f-point, where the a and vr bands are far from the 
vacuum energy [^. 

SV{z) describes the effect of higher-order image charges which arise at 0 < z < dg 
due to infinite multiple reflections of the primary image charges at the metal surface 
and the graphene layer. The latter can be regarded as metallic because of the efficient 
screening of lateral electric fields by the layer 29 . The classical image-potential 


between two metal surfaces with distance dg is given by 

11111 1 
4:Z 4:(dg — z) 4:dg 4:dg 4:(dg-\-Z) 4:{2dg — z) 




[Z = 


....( 2 ) 


The first two terms in ([^ represent the image-potential in atomic units caused by of 
the primary image-charges from each metallic surface. The following terms result from 
subsequent multiple reflections. By reordering, the latter can by written as 

1 1 


1 oo 2 
4 kdg kdg 


( 3 ) 


+ z {k l)dg — z 

This series converges to a net repulsive contribution. It increases the potential maximum 
in the center between the graphene layer and the metal surface (figure from 
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Figure 1. One-dimensional model-potential V{z) for grapliene/Ru(0001) at a 
graphene-metal distance of dg = 5.0 A. Black and red circles denote the position 
of Ru and C atomic layers, respectively. 


Vim(z = dg/2) = —l/dg to Vim{,z = dg/2) = 


the self-energy of a charged particle located in the gap between two metals 31 


ln(2)/(ig, which is the classic value of 

For 


decreasing distance dg, the potential maximum is lowered. However, even at dg = 2.2 A, 
as found in L areas of g/Ru(0001), ldm(dg/2) = —4.54 eV remains above the Fermi level. 
A simplihed description of the potential between metal and graphene by an effective 
quantum well with a large depth of 4.8 eV below the Fermi level, as proposed by Zhang 
, thus results in unrealistically large binding energies of the image-potential 


et al. p 

states. This is the main origin of the discrepancy in the identihcation of the hrst held- 


emission resonance observed in the STS spectra 19 -22 


The correction Vq,{z) considers the charge transfer between metal and graphene. It 
is approximated by a linear interpolation of the work function difference A<F between 
the clean metal surface and the combined system: 


14(z) = Ad) 1^1 - 0 < z < dg. (4) 

On the basis of the model potential V{z), which is depicted for the example 
g/Ru(0001) at a graphene separation dg = 5.0 A in hgure the wave function of 
the image-potential states and the corresponding binding energies with respect to 
the vacuum level F^vac have been obtained by solving the one-dimensional Schrodinger 
equation numerically by Numerov’s method. By extending the model potential to more 
than 100 A into the metal bulk and to more than 1000 A into the vacuum, we made 
sure that the results do not depend on the system’s extension. Even if we halve these 
intervals, the binding energies of the image-potential states with quantum numbers 
77, < 4 change by much less than 0.1 meV. 
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Figure 2. (a) Wave functions ^'(z) of the first two members of the symmetric (+) and 
antisymmetric (—) image-potential states of freestanding graphene, respectively, (b)- 
(f) Probability densities |'I'(z)p of the image-potential states at the graphene/Ru(0001) 
interface for graphene-Ru separations dg of 20 A (b), 10 A (c), 3.7 A (d), H areas) 
and 2.2 A (e), L areas) and at the clean Ru(OOOl) surface (f). The positions of the 
uppermost Ru atomic layer and the graphene layer are shown by vertical black dashed 
and red solid lines, respectively. The solid black line depicts the one-dimensional 
model-potential. The projected Ru bulk band structure (gray shaded areas) has been 
extended up to the image-plane of the metal {z = Zim)- 


3. Results and Discussion 

Table summarizes the calculated and experimental binding energies with respect 
to of the hrst two image-potential states for graphene on Ru(OOOl), Ir(lll), and 
Ni(lll). It shows that our model achieves a very good quantitative agreement with 
the experimentally observed binding energies which have been obtained by 2PPE for 
g/R(0001) and g/Ir(lll) [^. In particular, the distinct binding energies of the 
{n = 1) state in the L and H areas of g/Ru(0001) can be reproduced. Together with 
the predicted values for g/Ni(lll), this data suggests that the binding energies of the 
image-potential states are strongly correlated with the graphene-metal separation dg. 
In order to emphasize this correlation, we will discuss in this section the systematics 
of the formation of image-potential states at different graphene-metal separations for 
the example of the strongly corrugated g/Ru(0001) because its L and H areas exhibit 
distinct graphene-metal separations, that are typical for strongly and weakly interacting 
graphene/metal systems, respectively [I^ . 

The L and H areas of g/Ru(0001) form a hexagonal moire-superlattice with a lateral 
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Table 1. Calculated binding energies (eV) with respect to the vacuum level of 
the first two image-potential states (u = 1) and (n = 2) in the L and H areas of 
graphene/Ru(0001) and on almost flat graphene/Ir(lll) and graphene/Ni(lll). The 
results of the one-dimensional model-potential are compared to the experimental values 
obtained by 2PPE 
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23 


dg (A) denotes the respective graphene-metal separations. 



dg 

ID-Potential 

Experiment 

El 

£’2 

El 

£2 

S/Ru(0001) H 

3.7 

0.95 

0.24 

0.80 ±0.04 

0.18 ±0.06 

S/Ru(0001) L 1^ 

2.2 

0.47 

0.16 

0.41 ± 0.04 

0.18 ±0.06 

g/Ir(lll) g 

3.4 

0.88 

0.23 

0.83 ±0.02 

0.19 ±0.02 

g/Ni(lll) 

2.1 

0.41 

0.15 




32 


33 


Within our one- 


edge length of approximately 30 A of its 12.5 x 12.5-subunits 
dimensional model we can treat both areas separately because the lateral conhnement 
due to the hnite size of each area will reduce the binding energy of the image-potential 
states only slightly, as can be estimated from the area size of a half subunit A ^ 420 
to h^7r^/2mefrv4 43 meV, where m^s = 2.1 me with respect to the free electron mass 
ll. Since charge-transfer between the metal and graphene leads to a doping of 


Trip 


graphene 29,34 , we model the work function of g/Ru(0001) as an empirical function 
of the graphene-Ru separation dg by ‘hg/Ru(dg) = <hg — ^/dg, which makes it possible to 
vary dg systematically for the model calculations. 

Setting .^ = 1 eV A matches the work function $h = 4.24 eV in the H areas 


23 


as 


well as the 0.24 eV lower value in the L areas 29 . For large separation, it converges to 


<hg = ^(dg = oo) = 4.48 eV 35 of freestanding graphene. 

Figure shows the calculated probability densities of the image-potential states for 
different separations dg between the graphene layer and the Ru(OOOl) surface as well as 
for the separate systems. The corresponding binding energies are depicted in hgure 
For freestanding graphene, the two series with even (-b) and odd (—) symmetry can 
clearly be assigned from the amplitude of their wave functions at z = dg where the odd 
states have a node, whereas the even states have an antinode. The presence of the metal 
substrate breaks the mirror symmetry of the freestanding graphene and the parity of the 
image-potential is, in principle, no longer a good quantum number. At a graphene-metal 
separation as large as dg = 20 A, however, the influence of the substrate on the lowest 
members (n = !■*■, 1“) of the image-potential states is still weak due to the rather small 
spatial extent of their wave functions perpendicular to the surface. Therefore, their 
probability densities almost resemble those of freestanding graphene. We point out, 
that even for large dg, only image-potential states that originate from graphene and no 
states that are bound in the potential well formed by the image-potential of the metal 
and its surface barrier are found. This is related to the large difference in work function 
between graphene and Ru(OOOl) (<Fru = 5.51 eV 1^), which results in an additional 
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repulsion from the Ru surface. This contributes to a reduction of the binding energies 
with respect to those of freestanding graphene. The latter are depicted by dashed lines 
in figure]^ The same argument applies for the Ir(lll) (<f)ir = 5.79 eV [^) and the 
Ni(lll) (<hNi = 5.27 eV [^) surface. Beside the work function difference, the repulsion 
due to the surface barrier of the metal contributes additionally to a reduction of the 
binding energies with respect to freestanding graphene. At dg = 20 A, this applies 
particularly to the higher members (n = 2^, 3^, ...) because of their more extended 
wave functions, which are already considerably perturbed by the repulsion due to the 
surface barrier. This also results in an asymmetry and phase shift of the wave functions 
with respect to those of freestanding graphene. The wave function of the state labeled 
as n = 2“ at dg = 20 A, for example, has now almost an antinode at ^ = dg and more 
closely resembles the wave function of the {n = 3^)-state of freestanding graphene as 
denoted by the corresponding symbol in figure 

At intermediate separations, the properties of the image-potential states start to be 
governed by the image-potentials of both surfaces, and the attractive image-potential 
of the metal partly compensates for the repulsion due to the surface barrier and the 
work function difference. The attraction of the metal leads to a maximum of the 
binding energy of each state at specific separations as can be seen in figure This 
is also reflected by the form of the wave functions shown in figure |2 
for example, the symmetric (n = l+)-state is still almost unaffectec 
extended wave function of the (n = l“)-state is slightly attracted towards the metal 
surface which results in an increase of the probability density in between the metal 
surface and the graphene layer compared to larger separations. For comparison we have 
also made simplified calculations (not shown) where we described the metal as a simple 
barrier without an image-potential similar to the model used in 26 . In this case, we 


At dg = 10 A, 

while the more 


do not observe this attraction and all image-potential states are always repelled when 
the graphene layer approaches the metal surface. 

At small graphene-metal separations {dg < 5 A), the repulsion due to the surface 
barrier of the metal dominates and also the phase of the (n = 1“) wave function shifts 
and the latter starts to resemble that of the higher-lying (n = 2’'')-state of freestanding 
graphene. In particular at a separation of dg = 3.7 A, which corresponds to the H 
areas of g/Ru(0001), the vacuum part {z > dg) of all wave functions almost resembles 
the corresponding part of the symmetric series of freestanding graphene. This can be 
regarded as a transition from the double series of image-potential states of freestanding 
graphene to one single series, omitting the differentiation with regard to the parity. 
Since only the former {n = l’'')-state maintains a certain similarity to the corresponding 
state on freestanding graphene, we denote this fact by writing -|- in parentheses in the H 
areas. At a separation of dg = 2.2 A, which corresponds to the L areas of g/Ru(0001), 
all wave functions are further repelled from the metal surface. Even if the phase of the 
wave functions now correspond to the phase of the odd series of states on free standing 
graphene, the mirror symmetry has completely vanished, which is indicated by omitting 
the parity for the whole series. Now, all wave functions almost resemble the form of 
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Figure 3. Work function and calculated binding energies with respect to the vacuum 
level of the first four members of image-potential states at graphene/Ru as a function of 
the separation dg between the graphene layer and the Ru(OOOl) surface. Closed (open) 
symbols denote states with a symmetry at z = dg which is close to the symmetry of 
the even (odd) image-potential states of free-standing graphene. The binding energies 
of the first two even {n = 1+, 2+) and odd (n = 1“, 2~) image-potential states of 
freestanding graphene are depicted by solid and dashed horizontal lines, respectively. 
The vertical lines indicate the graphene-Ru separations of the L and H areas. The 
inset shows the experimental 2PPE spectrum adapted from 23 which averages over a 
large number of L and H areas. 


those on the clean Ru surface with one layer of atoms atop, as can be seen by comparing 
hgurej^e) and (f). In comparison to the clean Ru surface, however, the spatial extent of 
the wave functions in the vacuum is larger, which indicates a stronger decoupling from 
the surface. In the H areas, on the other hand, a considerable amount of the {n = 1*^+^) 
wave function still hts in between the graphene and the Ru surface and is therefore more 
strongly coupled to Ru as compared to the L areas. This state is closely related to the 
interlayer state of graphite [^. 

As depicted by the different symbols in hgure the change of the wave function’s 
character upon approaching the graphene layer to the metal surface is also reflected by 
the change of the binding energy with respect to the vacuum level. Whenever the binding 
energy of the image-potential states of the entire system becomes smaller compared to 
the next energetic adjacent state of freestanding graphene, the local symmetry at z = dg 
is closer to the latter. At dg = 20 A, for example, the binding energies of the hrst, second, 
and third state of the entire system are larger compared to the second, third, and fourth 
state of freestanding graphene. Consequently, their local symmetry at z = dg is still 
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close to that of the first (n = 1+), second (n = 1“), and third (n = 2+) state of 
freestanding graphene. The energy of the fonrth state of the entire system, however, is 
already smaller compared to the fifth state (n = S’*") of freestanding graphene. Its local 
symmetry at z = dg more closely resembles the symmetry of that state and not of the 
fourth state (n = 2~) of freestanding graphene. The energy of this state is therefore 
drawn as a closed symbol which denotes states of even (local) symmetry. For decreasing 
graphene-metal separation, the binding energies of the image-potential states increase 
slightly at first due to the attractive image-potential of the metal. At separations found 
in real graphene/metal systems (dg < 5 A), however, the repulsion due to the surface 
barrier of the metal dominates and results in a rapid decrease of the binding energies 
upon further approach. The sequence of the relative energy spacings is, however, rather 
insensitive towards dg and resembles a scaled Rydberg series with decreasing energy 
separation for increasing quantum number for all dg. 

The absolute change of the binding energy as a function of dg is most pronounced 
for the first image-potential state. It decreases from more than 1.4 eV at large dg to 
0.95 eV at the separation found in the H areas of g/Ru(0001) and finally to 0.47 eV at 
the separation found in the L areas. The difference between the binding energies with 
respect to i^vac in the H and L areas of 0.48 eV is twice as large as the difference between 
the local work functions of 0.24 eV 29 . The first image-potential state in the H areas 


has therefore a smaller energy with respect to the Fermi level when compared to the 
the corresponding state in the L areas. This confirms the assignment of the peaks in the 
2PPE spectrum of which is reproduced in the inset of figure This assignment i 


IS 


also in agreement with density functional theory (DFT) calculations by Borca et al. 19 


which shows that the results of our model calculations are robust also for small distances 
where DFT can explicitly describe the chemical interaction between graphene and the 
metal. For large distances and higher quantum numbers the evolution of the binding 
energies with distance depicted in 19 considerably differs from our results which we 
attribute to the fact that DFT cannot properly consider the long-range image-force. 

For the second image-potential state, the calculated difference between binding 
energies with respect to Evac in the H and L areas is in the order of the difference in the 
local work function, which explains that these two states could not be resolved separately 
in the 2PPE spectra. The red horizontal lines in figure depict the experimental results 
from 23 under consideration of the local work functions in the H and L areas, and 


help visualize the very good agreement between experiment and model calculation. The 
small deviation for the binding energy of the lowest image-potential state energy in the 
H areas is reduced if one considers the additional reduction of the calculated binding 
energy due to the lateral localization as discussed above. The localization of this state 
should be stronger in the H areas as deduced from the effective mass which has been 


found to be twice as large as the free electron mass 23 


The results of our model calculation make it also possible to qualitatively interpret 
the experimentally observed decay dynamics of electrons in these states, as has been 
investigated by time-resolved 2PPE 23 . As shown in figure [^e), the probability 
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densities of all image-potential states in the L areas mostly resemble those of the 
clean Rn(OOOl) snrface bnt with some degree of deconpling, i.e. the center of gravity 
of the probability density is located further away from the metal into the vacuum. 
At hrst glance, this similarity to the clean Ru surface seems to be surprising. The 
graphene-metal separation in the L areas is, however, almost identical to the Ru-layer 
spacing (^ru = 2.14 A 37 ). The Ru surface covered by one metallic-like graphene 
overlayer therefore has similar properties with respect to the image-potential states as 
the extension of the Ru crystal by one more layer of Ru atoms. At the clean Ru(OOOl) 
surface, the lifetimes of electrons in the hrst two image-potential states of Ti = 11 fs 
and T 2 = 57 fs 38,39 are comparably short because the decay by inelastic scattering 
with bulk electrons is very efficient due to the high density of states of the Ru 4d bands 
just below the Fermi level 39 . Thus, the slightly enhanced experimental lifetimes of 
Ti- = 16 ± 5 fs and T 2 - = 85 ± 13 fs 23 are a direct indication for a certain degree 


of decoupling in the L areas. In the H areas, a considerable part of the probability 
density of the {n = l^+^j-state hts below the graphene layer close to the Ru surface. 
This enhances the electronic coupling to the Ru bulk states as is rehected by the shorter 
lifetime of Tp+) = 11 ±8 fs 23 , which is identical to the value found on clean Ru(OOOl). 

By adapting Rm(^), our model can be easily applied to the description of graphene 
layers on other metal substrates. With Ir(lll) and Ni(lll) we chose two examples 
that represent limiting cases of weak and strong interaction with the graphene layer, 
respectively. It has been shown that g/Ir(lll) grows with a rather large and 
homogeneous separation of dg = 3.4 A 


40 


which is very similar to the separation 
found in the H areas of g/Ru(0001). Graphene and Ni(lll) have almost no lattice 
mismatch. This results in the growth of a flat graphene layer on Ni(lll), despite the 
strong interaction 16 . The resulting separation of about dg = 2.1 A is very similar 
to that found in the L areas of g/Ru(0001). Binding energies and lifetimes of the 
image-potential states on g/Ir(lll) have been measured by 2PPE 11 and can be 
directly compared with the results of our model calculations. No experimental data 
on g/Ni(lll) are available, but 2PPE results on the binding energies for the clean 
Ni(lll) surface 36 make it possible to properly adjust the parameter of the metal 
potential and to predict the binding energies for the graphene covered surface. The 
ferromagnetic coupling of Ni results in a spin-split surface-projected band structure 
which is reflected by an exchange-splitting of the image-potential states. The observed 
splitting of 14 ± 3 meV is, however, small compared to the linewidth and can only be 
observed with spin-resolved detection 36 . We therefore neglect the exchange-splitting 


here and use a spin-integrated surface-projected band gap extracted from [4^ for our 
model calculations. Because the size and position of the surface-projected bulk band 
gap do not differ much between Ru(OOOl), Ir(lll), and Ni(lll), the calculated binding 
energies depend most sensitively on the graphene-metal separation dg. Consequently, 
the binding energies for g/Ir(lll), as well as for g/Ni(lll), are very similar to the results 
for g/Ru(0001) at the corresponding separations. As shown in table the results of 
our model calculation for the binding energies on g/Ir(lll) are in excellent agreement 
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with the experimental data. For g/Ni(lll) we predict much smaller binding energies 
compared to g/Ir(lll). This is due to the stronger repulsion of the image-potential 
states at the smaller separations. The binding energies on g/Ni(lll) are on the other 
hand comparable to those on the L areas of g/Ru(lll) in correspondence to the similar 
separation dg. 

Finally, we would like to comment on the formation of image-potential states at 
the interface of graphene/SiC, even if our model is not directly applicable to this 
semiconducting substrate. Although graphene might possess only a weak interaction 
to SiC, our results generally show that the shape of the image-potential states is very 
sensitive to the substrate induced symmetry break even for large separations due to 
their large spatial extent. In another 2PPE experiment on graphene/SiC(0001), Shearer 
et al. 42 indeed observed a single series of sharp symmetric peaks of image-potential 


states in agreement with data shown in [^. Therefore, it seems rather unlikely that the 
experimentally found splitting of the first image-potential state on graphene/SiC [8,10 


can be interpreted as a remnant of the mirror symmetry of freestanding graphene with 
a small energy separation between these two states. The latter is particularly hard 
to understand because the approach of a flat graphene layer towards a substrate with 
a surface barrier does not change the relative energy spacing of adjacent states, i.e. 
the energy spacing between the lowest two states is always larger than the subsequent 
spacing because the energy sequence of the Rydberg series is at all distances dominated 
by the long-range image-potential. 


4. Conclusion 

On the basis of an analytical, one-dimensional potential, we can quantitatively describe 
the formation of image-potential states at the interface between a graphene layer and a 
metal surface by means of model calculations. By systematic variation of the graphene- 
metal separation, we have shown how the double Rydberg-like series of even and odd 
image-potential states of freestanding graphene evolves towards a single series of the 
compound system when a flat graphene layer is located at a distance to the metal as 
found in real graphene/metal systems. This transition is driven by the repulsion from 
the metal substrate, which increasingly reduces the mirror symmetry of freestanding 
graphene with decreasing separation. Even for separations found in weakly interacting 
systems, we find no remnant of the distinction between states of even and odd symmetry. 
At these intermediate separations, the first image-potential state is partly trapped 
between the metal and the graphene layer. It attains properties of an interfacial 
state. At distances found in strongly interacting systems, on the other hand, the wave 
functions almost resemble those of image-potential states on a clean metal surface but 
with a higher degree of decoupling. The repulsion of the image-potential states from 
the metal surface is also responsible for a strong reduction of their binding energies 
with decreasing graphene-metal separation. This explains the distinct differences of the 
experimentally observed binding energies of the first {n = 1) image-potential states in 
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the H and L areas of strongly corrugated g/Ru(0001). According to the respective 
similar separations on weakly interacting g/Ir(lll) and strongly interacting g/Ni(lll), 
we hnd respective comparable binding energies which are in excellent agreement with 
available experimental data for g/Ir(lll). 
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Appendix A. Metal Potential 


For the metal part V^{z) of our model potential we used the well-established one¬ 
dimensional potential introduced by Chulkov et al. 
is piecewise dehned along the surface normal by: 


This parameterized potential 


Yrr^iz) = 


—A-io + Al cos ( ——^ 
\dm 


-A 20 -F A 2 cos(fiz) 
-A 3 exp[-a ;(;2 - Zi)] 
exp[-A(2; - Zira)] - 1 


4(z 


z<t) 

0 < z < zi 

Zi < Z < Ziin 
Zin. < Z 


(A.l) 


Here dm is the layer spacing and z = 0 corresponds to the position of the surface 
atom. By requiring that Vm{z) and its hrst derivative to be continuous at the matching 
points Zi and Zim, only four of the ten parameters Ai, Aio, A 2 , ^ 20 ; A 3 , a, /3, A, Zi, 
Zim are the independent of each other. It is feasible to choose the offset Ai and the 
amplitude Aio in order to reproduce the energetic position and width of the surface- 
projected band gap. For Ru(OOOl), Ir(lll) and Ni(lll) (spin averaged) these have been 


extracted from 43, 44] , [IT] and |4l] , respectively. The parameters A 2 and /3 have been 


used to reproduce the experimental binding energies En of the image-potential states. 
For Ru(OOOl) we have htted these parameters to the experimental binding energies 
El = 0.66 eV and E 2 = 0.19 eV reported by Gahl et al. 38 . For Ni(lll) we have 
used data reported by Andres et al. 1^. Due to the lack of experimental data on the 


image-potential states of clean Ir(lll), we have estimated the binding energies by using 
the Rydberg formula En = 0.85 eV/(n -|- a)^ where a is the quantum defect. For a sp- 
inverted band gap, the latter can be determined from the position of the vacuum energy 
Evslc within the projected band gap [^. This results in Ei = 0.64 eV and E 2 = 0.18 eV. 

The parameter sets used for modelling Ru(OOOl), Ir(lll), and Ni(lll) are collected 
in table I 
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Table Al. Parameters of the metal part Vmiz) of the model-potentials used for 
Ru(OOOl), Ir(lll), and Ni(lll). The work function $, the potential offset Aio and 
amplitude Ai are given in eV, the layer spacing in A and the parameter /3 in A“^. 



dm 


$ 


dio 

Al 

^2 

P 

Ru(OOOl) 

2.138 

37 

5.51 

23 

7.436 

9.400 

12.70 

4.3464 

Ir(lll) 

2.217 

45 

5.79 

11 

9.511 

6.188 

6.50 

4.6068 

Ni(lll) 

2.035 

46 

5.27 

36 

11.259 

6.503 

4.91 

2.3070 



Position z (A) 


Figure Bl. One-dimensional model-potential of freestanding graphene. The vertical 
dashed lines delimit the different definition intervals. The red dashed line represents the 
asymptotical form of the classical potential due to the image-charges of both surfaces 
(inset). Gray circles depict the ”LDA-|-image tail” potential which has been extracted 
from figure 1 of [^. 


Appendix B. Graphene Potential 


Our analytical approximation of the ”LDA-|-image tail” hybrid potential of Silkin et 
al. for the description of freestanding graphene is inspired by the dehnition of the 
metal potential with omission of the bnlk part: 



"S 

1 

0 

-b Gi cos(/3i 2 ;) 

0 < 

\z\ 

< 

kll 



1 

to 

0 

- G2COs[I32{\z\ - |2:i|)] 

kil 

< 

z\ 

VI 


UP) = 

-G 3 exp[-a(| 2 ;| - \z2\)] 

1 ^ 2 ! 

< 

z\ 

^ 1 Zira 1 

(B.l) 


exp 

-A(|z| - l^iml)] - 1 


< 






4(1^1 - |2:im|) 

1 Avu 

\z 




As shown in hgnre Bl, Vg{z) is symmetric in the 2 ;-direction with respect to the 
graphene layer aX z = 0 and converges to the classical image-potential for large distances 
1^1. Again, only fonr of the twelve parameters Zi, Z 2 , Gi, Giq, G 2 , G 20 , G^, /3i, 
/? 2 , and A are independent if we reqnire Vg{z) to be continnonsly differentiable at 
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Table Bl. Binding energies (eV) calculated with the one-dimensional model 
potential and the ”LDA-|-image tail” hybrid potential [^, respectively. 



Ei+ 

El- 

± 2 + 

± 2 - 

ID-Potential 

1.46 

0.60 

0.28 

0.18 

LDA-|-image tail 

1.47 

0.72 

0.25 

0.19 


the matching points Zi, Z 2 , and z-^^. We choose the potential offset Gio = Vq, 

the amplitude Gi, and the inverse widths fli and /32 to fit Vg{z) to the ”LDA-|-image 
tail” hybrid potential for a matching point between LDA potential and image-tail of 
Zq = 3 a.u. The best fit is shown as solid line in figure |B1[ It has been achieved with 


Gi = 3.336 eV, Gio = 29.3465, /3i = 8.7266 and = 14.8353. 

The binding energies Ef; of the symmetric and antisymmetric image-potential states 
calculated with these parameters are listed in table ^ They agree well with those 
calculated using the ”LDA-|-image tail” hybrid potential for zq = 3 a.u. The two 
equivalent positions of the image-planes of z^^ = ±0.99 A, which result from the chosen 
parameters, are almost coincident with the spatial extent of the polarizable conjugated 
TT-system [^. 

The screening of electric fields by the metal as well as the graphene layer is taken into 
account by cutting the corresponding image tail of \fa{z) and Vg{z) at the respective 
opposite surface {y^{z) = 0 for 2 ; > dg and Vg{z) = 0 for 2 ; < 0). Together with 
the corrections SV(z) and 14-(z), this results in a discontinuity of the total model 
potential V(z) at z = 0 and z = dg. The discontinuity is compensated by increasing the 
amplitudes of the respective cosine oscillation of the metal potential Idn (A.l) and the 
graphene potential 14 (B.l). 
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